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ABSTRACT 

We have used XMM-Newton to observe six Lyman Break Analogs (LBAs): 
members of the rare population of local galaxies that have properties that are 
very similar to distant Lyman Break Galaxies. Our six targets were specifically 
selected because they have optical emission-line properties that are intermediate 
between starbursts and Type 2 (obscured) AGN. Our new X-ray data provide 
an important diagnostic of the presence of an AGN. We find X-ray luminosities 
of order 10^^ erg s~^ and ratios of X-ray to far-IR luminosities that are higher 
than values in pure starburst galaxies by factors ranging from ~ 3 to 30. This 
strongly suggests the presence of an AGN in at least some of the galaxies. The 
ratios of the luminosities of the hard (2-10 keV) X-ray to [O III]A5007 emission- 
line are low by about an order-of-magnitude compared to Type 1 AGN, but are 
consistent with the broad range seen in Type 2 AGN. Either the AGN hard 
X-rays are significantly obscured or the [O III] emission is dominated by the 
starburst. We searched for an iron emission line at ~ 6.4 keV, which is a key 
feature of obscured AGN, but only detected emission at the ~ 2a level. Finally, 
we find that the ratios of the mid- infrared (24/im) continuum to [O III]A5007 
luminosities in these LBAs are higher than the values for Type 2 AGN by an 
average of 0.8 dex. Combining all these clues, we conclude that an AGN is likely 
to be present, but that the bolometric luminosity is produced primarily by an 
intense starburst. If these black holes are radiating at the Eddington limit, their 
masses would lie in the range of 10^ to 10^ Mq. These objects may offer ideal 
local laboratories to investigate the processes by which black holes grew in the 
early universe. 



Subject headings: galaxies: starburst — quasars: general — X-rays:galaxies 
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Introduction 



One of the major unsolved problems in astrophysics today is the nature of the connection 
between the formation and evolution of galaxies and supermassive black holes (SMBHs). 
Current models invoke a close interplay ( "feedback" ) between SMBHs and their host galaxies 
in order to arrest the growth of the mos t massive galaxies, which are otherwise over-predic ted 
by models of structure formation (e.g., ISpringel fc Hernquist II2003I : ICroton et al. 1120061 ). 



Observations of high-redshift galaxies paint a somew hat confusing pic ture of the relation- 
ship between star formation and black hole growth (e.g.. iHeckman 1120091 ). The best-studied 
population of high-redshift star-forming galaxies are the Lyman-Break Galaxies (LBGs) and 
closely related types of UV-selected galax ies. They repres ent a major phase in the early 
stages of galaxy formation and evolution (iGiavalisco II2002I ). They are a primary tracer of 
high redshift star formation, as they are easily found aX z = 2 — 6 in deep pencil beam optical 



surveys (e.g.. iSteidel et al. Ill996l : IShaplev et al. II2003I : IGiavalisco et al. 1120041 ). The LBGs 
are believed to be the precursors of present day galaxies undergoing a phase of intense star 
formation and merging, and a large fraction of LBGs may merge to form early-typ e galaxies 
that are situated in present-day groups and clusters (e.g., lAdelberger et al. 1120051 ). 



Most of the LBGs do not harbor an Active Galactic Nucleus (AGN) luminous enough 
to be individually detected in the deepest X-ray sur veys. However , a subset (~3%) of LBG s 



with actively accreting SMBHs has been identified ( ISteidel et al. II2002I : iLaird et al. II2006I ). 



These AGNs typically have X-ray luminosities of 10^^ — 10^^ erg s~^ in the 2-10 keV energy 
range, and often show signs of significant obscuration. Detailed studies of what might trigger 
both the AGNs and the starbursts in these forming galaxies are difficult since their distances 
render them small and faint. 

The Galaxy Evolution Explorer (GALEX) mission is conducting surveys in the far- and 
near-UV (A ^ 1600 and 2300 A respectively) and has been used to identify a rare popula- 
tion of low-redshift {z ~0. 1-0. 3) galaxies that are rernarkably similar to high-redshift LBGs 



in ni o st physical properties ( jHeckman et al. I l2005al : iHoopes et al. 



2008 



2009 



20101). This population was first identified by 



2007 



Overzier et al. 



Heckman et al. 



(l2005a[ ) and 



Hoopes et al. " (EoOT ). based on cross-matching the GALEX and Sloan Digital Sky Survey 
data. These objects were defined by two criteria designed to match the typical LBGs: 1) a 
minimum luminosity in the far-UV (FUV) band Lpuy > 2 x lO^^L©; 2) an FUV effective 
surface brightness Ipuv > lO^LQkpc^^, where Ipuv is the mean FUV surface brightness 
interior to the SDSS u half-light radius 0. These local "Lyman break analogs" (LBA) are 



^While selection based on high FUV luminosity and surface brightness also selects Type 1 AGN and BL 
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similar to distant LBGs in size, surface brightness, metallicity, dust, star formation rate 
(SFR), stellar mass, and gas velocity dispersion. There are other local star forming galaxies, 
e.g., Ultra-Luminous Infrared Galaxies (ULIRGs) and Blue Compact Dwarfs (BCDs), which 
have some similarities with LBGs. However, ULIRGs are much more heavily obscured by 
dust, and BCDs are much smaller in mass and luminosity (star-formation rate) compared 
with LBGs. These differences make them poorer local analogs to LBGs. However, it is 
intriguing that a n actively accreting SMBH has recently been reported in a nearby dwarf 
star burst galaxy (IReines et al. II2OIII ). 



Interestingly, about 20% of the LBAs have optical emission-line spectra that are inter- 
m ediate between those of pure starbursts and those of Type 2 (obscured) AGN. As discussed 
by lOverzier et al. I (120091 ). the optical spectra alone do not unambiguously establish the pres- 
ence of a Type 2 AGN. They may instead be extreme starburst systems in which the physical 
conditions in the ionized g different from typical local starbursts. 



Overzier et al. I (120091 ) also presented imaging with the Hubble Space Telescope that 



showed that most of these possible starburst-AGN composites contain extremely massive 
(one to few-billion Mq) and compact (radius ~ 10^ pc) starbursts. These Dominant Compact 
Objects (DCOs) would appear to be ideal sites for the formation and/or rapid growth of 
a SMBH. It has been suggested that the seeds of present-day SMBHs may be present at 
the c enters of the massive, compact stellar clum ps seen in clumpy galaxies at high redshift 



[e.g. lElmegreen. Bournaud &: Elmegreen II2008I ). If these LBAs indeed host an AGN, they 
would offer the opportunity to study the relationship between the formation of galaxies 
and supermassive black holes in detail in objects that appear to be excellent local analogs 
to a major population of high-redshift galaxies. Unfortunately, the optical spectra alone 
do not unambiguously establish that these LBAs are indeed compos ite starburst /obscured - 
AGN systems (as opposed to starbursts with extreme properties) (lOverzier et al. II2009I ). 
Observations in the hard X-ray band can potentially offer us important new diagnostics of 
the presence of an obscured AGN. 

In this paper we report the results of such an investigation. We describe our data in 
§2] and present the spectral analysis in §31 The results are discussed in §11 and we come to 
some brief conclusions in §S1 



Lacs, the se objects have been el iminated from samples of LBAs (see lHoopes et al. 1120071 for details). Note 



also that Overzier et al. ( 20091 ) presented very high signal-to-noise optical spectra of LBAs showing that 
there is no spectroscopic signature of a Type 1 AGN. 
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2. Data Analysis 



2.1. Sample selection 

To select our sample f or X-ray observations we begin with the 31 LBAs with HST imag - 
ing (jOverzier et al. 1120091 ) and Spitzer (IRAC+ MIPS) photometry (jOverzier et al. Il201ll ). 



In Fig. [T]we show the standard "BPT" diagram (IBaldwin et al. Ill98ll ) used to classify the 



optical emission-line s pectra of galaxies, where the emission line ratios are from Table 2 in 



Overzier et al. 



(|2009[ ). We plot all the galaxies in the GALEX-SDSS cross- matched sample 

and then highlight the 31 LBAs in our sample (large symbols). We note that sevenE of the 
LBAs lie in the composite region in the "BPT" diagram, which is known to harbor "com- 
posite objects" , which ar e typically (but not nece s sarily) objects havin g both star-formation 
and a Type 2 AGN (see iKauffmann et al. II2003I : iKewley et al. II2006I ). Hereafter, we refer 
to this sample as "LBA-composites" . In this paper, a fiat A cold dark matter (ACDM) 
cosmology with h = 0.7, = 0.7 is assumed. 



2.2. The X-ray data 

AMM-Newton observations were performed for 6 of the 7 LBA-composites. SDSS 
J0054+1554 was given a "C" priority and not observed (marked as green asterisk in Fig. [T]). 
A summary of observations is shown in Table [H 

The reduction of the AMM-Newton data was done by XAssist (ver . . 97)1, which 
is a software for automatic analysis of X-ray data. XAssist runs the AMM-Newton Science 
Analysis System (SAS) packages automatically to filter the data, generate the light curves and 
extract the spectra of sources in the field, as well as compute the associated redistribution 
matrix files (RMF) and ancillary region files (ARF). The LBA-composites are weak in X-ray 
emission, and their optical positions were used to define the source regions in X-ray data 
reduction. XAssist determines the source extraction sizes for source by fitting an elliptical 
Gaussian model to the source image, which in the case of on-axis AMM observations typically 
results in regions of size 20-25". Background regions were chosen as annuli centered on the 
sources to extract the background spectra. Table H] lists the PN and MOS counts for each 
source. 



^We include SDSS J1434+0207 (red cross in Fig. [T]) in our "LBA-composite" sample, although it lies 
close to the lower left of the solid line in the BPT diagram and does not harbor a DCO. 

■^http: / / xassist.pha.jhu.edu/ zope/xassist 
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2.3. Ancillary data 



The Optical Monitor (OM) aboard on XMM-Newton allows us to make simultaneous 
observations in optical/UV and X-ray bands. We processed the OM data with the omichain 
package in SAS to obtain UV count rates and the corresponding magnitudes of each object, 
as well as their associated uncertainties. We obtained exposures in the U (344nm), UVWl 
(291nm), UVM2(231nm) and UVW2 (212nm) band filters, but only the UVWl was used 
for all six LBAs. We list the UV detections with each filter and the luminosity measured 
in UVWl in Table |2l The conversions from UV count rates to specific fluxes (in units of 
erg/s/cm^/A) were calculated based on the white dwarf standard star^. 

Table[3]lists the observed and extinction-corrected [O III]A5007 line luminosity (L[oiii],obs 
and -L[oiii],corr, respectively) for each LBA-composite. We obtain L[oiii],obs from the SDSS 
Data Release 7 MPA/JHU catalog. We derive i^[oiii],corr from -L[oiii],obs using the observed 
Ha/H/3 ratio, an intrinsic (Ha/H/3)n rati o of 2.9, and the standard extinction curve for 
galactic dust ( lOsterbrock fc Ferland Il2006[ ). In Table El we also list the mid-infrared ( MIR) 
and far-infrared (FIR) luminosity, Iv24^m and Lpm, respectively (lOverzier et al. 1120111 ). 



3. Spectral Analysis 
3.1. X-ray Spectra 



X-ray spectral fitting was performed with XSPEC (ver. 12). X-ray photons are collected 
by three detectors on AMM-Newton, i.e. PN, MOSl and M0S2. The data of three detectors 
can be fitted simultaneously in XSPEC by tying all parameters together except a multiplicative 
constant factor. We froze this factor at 1 for PN, and allowed it to vary between 0.8 to 1.2 
for MOSl and M0S2. Since our targets have limited counting statistics (see Table S]), we 
first fit their spectra using the C-statistic ( ICash Ill979l ). which is often used in this low counts 
regime. The spect ra were grouped to 1 count per bin, to improve the performance of XSPEC 
( iTeng et al. 1120051 ). The simplest spectral model often fit to AGN is an absorbed power-law 
model. The X-ray photons suffer from two sources of absorption. The first is due to the 
neutral hydrogen in our own galaxy, and its column density (Ah,g) is listed in Tabled! The 
other is the intrinsic obscuration from the gas at the redshift of the source (AH,in)- AH,in 
is a free parameter in the spectral fitting as well as the photon index F of the power law. 
However, due to the low signal-to-noise ratio (the number of photons detected by the PN 



^http: / /xnim.esa.int / sas / current / watchout / 
Evergreen_tips_and_tricks/uvflux.shtml 
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and MOS detectors are listed in Table H]), we were not able to constrain both A/^n^n and F 
simultaneously in several cases. Therefore for sources with low S/N, we alternately froze 
iVn^in and F and let the other parameter vary. The same methodology was applied when 
adding an additional thermal component into the power-law continuum. 



3.1.1. Absorbed power-law model 

As mentioned above we fitted the spectrum in two ways to derive the power-law photon 
index and the intrinsic absorption: First, we freed the power- law photon index and assumed 
no absorption at the redshift of the source (A^H,in = cm~^) and that the continuum power- 
law spectra is only absorbed due to the Galactic column density (listed in Table [1]). Again 
in many cases there was not sufficient signal for the absorption and slope to both be free 
parameters. We list the photon indices derived from this simple absorbed power-law model 
in column 3 in Table IH as well as the corresponding C-statistic in column 4. Second, the 
spectra were fitted with the power-law model this time with the intrinsic absorption as a 
free parameter while the photon index was fixed as F = 1.9, which is a mean value for 



X-ray detected LBGs (iLaird et al. Il2006[ ) and is also a typical spectrum for an unobscured 
AGN. The values of A^H,in are shown in column 5 in Table ID The errors of photon indices 
in column 3 and the upper limits of intrinsic absorptions in column 5 were calculated at a 
confidence level of 90% (i.e., single interesting parameter AC=2.7). We also list in Table 
m their derived luminosities in soft X-ray (0.5-2.0 keV) and hard X-ray (2-10 keV) bands 
(column 7 and 8, respectively), and the ratio of hard X-ray and [O III] luminosities (column 
9). The associated uncertainties for luminosity correspond to the 90% confidence range of 
normalization of the continuum power-law component in spectral fit. See Fig. |2] for the 
spectral fits of each source with a fixed photon index F = 1.9. 

We also fitted the X-ray spectra of our targets using the x^-statistic to see how different 
the fits are compared with those by C-statistic. The fits were performed in the same way as 
those using the C-statistic above. The spectra were grouped in 10 counts per bin for both 
FN and MOS detectors, and the spectra of MOSl and M0S2 were combined in the fits. 
Table [5] lists the photon index and intrinsic column density of each source derived while the 
other one is frozen in the spectral fitting. SDSS J0802+3915 is excluded due to insufficient 
bins for fitting after grouping into 10 counts per bin. Also, spectral plots are shown in Fig. 
El which are fits with photon indices fixed at 1.9. Fitting the spectra using does not make 
any improvements or significant difference in estimating the absorption and photon index 
compared to the C-statistic method. Thus, we use the X-ray luminosity derived from the 
C-statistic fit in the following discussion. 
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Additionally, as SDSS J0808+3948, SDSS J1434+0207 and SDSS J0922+4509 have 
relatively more photons detected, we also fitted their spectra by thawing both the parameters 
of photon index and intrinsic absorption (using both the C and statistics). Table E] shows 
the values of these parameters and their associated errors or upper limits. Since we have two 
interesting parameters in the fits, the errors at the 90 per cent confidence level are calculated 
according to Ax^ or AC=4.6. 



3.1.2. Power-law plus thermal model 



The LBAs are the sites of strong starbursts, so we would expect to see thermal X-ray 



emission from the hot gas that is a characteristic features of such systems (e.g.. lGrimes et al. 



20051 . 12OO61 120071 ). We therefore fit models in which we added a thermal component to the 



absorbed power-law model aforementioned. We used the therni a l plas ma model apec in 
XSPEC instead of the older model, vmekal, used in iGrimes et al. I (120071 ). The apec model 
has three free parameters: the plasma temperature, the metal abundance (in solar units) and 
the normalization. Due to the low number of photon counts, we fixed the abundance at the 
solar value. The parameter of temperature and normalization were free in fitting. However, 
due to insufficient photon counts, we fixed the temperature at 0.5 keV for SDSS J0808-I-3948, 
SDSS J0922-t-4509 and SDSS J08 02+3915 which is a typical value for starburst galaxies (e.g.. 



Grimes et al. II2005I . l2006l . 120071 ). Here we kept the photon index fixed at 1.9 and allowed 



the intrinsic absorption at the redshift of the LBA to vary. The results are shown in Table 
[TJ which lists the temperature of plasma, the column density, and the derived luminosities of 
both the thermal and power-law component, as well as the C-stat for comparison with the 
single power-law fitting. The errors for parameters were calculated according to AC = 4.6, 
and the uncertainty of luminosity was proportional to the corresponding error of normaliza- 
tion of each component. For the three LBAs where the temperature was permitted to vary 
the temperature ranged between 0.3 and 0.6 keV. The thermal component luminosities were 
on the order of 10^^ erg s-\ with 10^° erg s"^ for SDSS J0808+3948 and about 10^^ g^g s"^ 
for SDSS J0802+3915. The power-law component luminosity in the hard X-ray range (2-10 
keV) is within a factor of 2 of those derived from the single power-law model spectral fit 
mentioned above. This power-law plus thermal model fit makes no statistically significant 
difference in the derived column density compared to that of the single power-law model. 

Since the signal-to-noise of the XMM-Newton spectra is low, we attempted to constrain 
any thermal component by simultaneously fitting the spectra of all LBAs together assuming 
that they have approximately the same hot gas temperature, i.e., we tied the temperature 
parameters together between the fits but let the apec normalization vary for each LBA. We 
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also fixed the power-law index at 1.9, and tied the normalization of power-law component in 
our fit (although note that overall normalization for each LBA was allowed to be free). The 
temperature given by this joint fit is 0.56 keV, and its 90% confidence level is 0.43—0.69 keV. 
We list the intrinsic column density, the thermal luminosity and the soft and hard X-ray 
luminosity contributed by the power-law component in Table [HI We see that the column 
densities and thermal luminosities given by this joint fit are consistent with those in the 
individual fits, but with improved constraints. 

In Fig. m we compared the two compositions of soft X-ray luminosity, i.e., the thermal 
and power-law luminosities in 0.5-2.0 keV, where the dashed line indicates where both values 
equate. This comparison was done separately for both the cases in which the spectra were 
fit individually and those that were fit simultaneously. These are shown in the left and right 
panel of Fig. H] respectively. The results derived from the individual spectra fits imply that 
these two mechanisms play a nearly equal role in soft X-ray emission, except in J0808+3948 
where the power-law soft X-ray luminosity is about one dex higher than that of thermal 
emission. For the joint spectral fit, the luminosities from the two emission mechanisms are 
comparable, within the range of 1 dex. The obscuration of SDSS J0214-I-1259 is Compton- 
thick in the joint fit, and the power-law luminosity at 0.5-2 keV is negligible (not shown). 



3.2. Star formation rates 



We list the UV luminosity for each LBA measured from the OM observation in Table [21 
We de rive the SFR from our UV luminosities using the formula from [Rosa-Gonzalez et al. 
(|2002[ ). which is based on evolutionary synthesis models of starbursts by assuming a Salpeter 
initial mass function (IMF). For consistency with [Overzier et al. I (120091 ). we divide the Rosa- 



5 to reflect the use of a Kroupa rather than a Salpeter 



Gonzalez et al. relati on by a factor of 

IMF (as described in [Calzetti et al. [ (|2009[ )). The resulting formula is as follows: 



SFR(UV) = 9.3 X 10-2^L^(erg s'^Hz-^) Mq/ji. 



In Table [91 we compare our SFR derived from the UV luminosity with that derived 
from a combination of the 24/im and Ha luminosities (see [Overzier et al. [[20091 for details). 
Due to the unce r tain i nfrared SEDs, no /^-correction was applied to the 24/im luminosity in 
Overzier et al. [ ( |2009[ ). and the actual SFRs might be somewhat higher than their derived 
values. The UV-derived SFR is lower than that derived from the IR continuum plus optical 
emission lines. This is consistent with the fact th at the UV photons su ffer significant dust 



attenuation and are re-emitted in the IR (see also [Overzier et al. [[201l[ ). 



- 9 - 



We note that any AGN contribution to the UV continuum is neghgible in these ob- 
j ects. This is demonstr ated most directly by HST COS FUV spectra of two of our targets 
faeckman et al. Iboilj ). 



4. Discussion 

As we described in the introduction, the optical emission-line spectra of our sample of 
LBAs suggest that these may be composite objects co nsisting of an intense s tarburst and an 



obscured AGN. However, this evidence is ambiguous. lOverzier et al. I (|2009[ ) have discussed 
this in detail and suggest as an alternative that the optical spectra of these galaxies could 
reflect extreme conditions in the ionized gas related to the presence of a starburst with 
extreme properties. 

The X-ray data we have presented provide important new insight into the nature of 
these objects. We begin by comparing the amount of detected hard and soft X-ray emission 
in the LBAs compared to pure starbursts with the same star-formation rates (as traced by 
the far-IR luminosity). We will then compare the LBAs to typical Type 1 and Type 2 AGN 
in terms of the ratios of the X-ray to [O III] and mid-IR to [O III] luminosities. 



4.1. X-ray properties compared to starbursts 



The far- infrared (FIR) luminosity is co nsidered to be a reliab le indicator of the global 



SFR of star forming galaxies. In particular, iRanalli et al. I ( 120031 ) found that the 2-10 keV 



X-ray luminosity has a tight linear correlation with the FIR luminosity by investigating a 
sample of local star forming galaxies. We show our LBAs in the L2-10 keV — -^fir diagram 
(see Fig. [7]) with the Ranalli relation, which is indicated by the solid line with two dashed 
lines indicating a factor of 2 above and below. The hard X-ray luminosities of the LBAs 
are taken from Table HJ which are derived from the absorbed power-law model fit with the 
photon index fixed at 1.9. For compar i son, w e al so show the sample of Luminous Infrared 
Galaxies (LIRG) from llwasawa et al. I (120091 ) and iLehmer et al. I (120101 ). where Iwasawa et 
al's sample was divided into two groups of AGN and "Hard X-ray Quiet" galaxies (HXQ) 
indicated by diamond and asterisk symbols respectively, and Lehmer et al's sample was 
indicated by triangle symbol^. We can see that our LBAs are about 0.5 to 1 dex above the 



iLehmer et al. I (|2010l ) suggest that the X-ray emission from high-mass X-ray binaries may be suppressed 
by heavy extinction in some of the HXQ, LIRGs and ULIRGs 
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Ranalli et al. relation, and even more displaced from the non-AGN LIRGs. This implies that 
a normal star forming process would only account for a minority of the hard X-ray emission 
of these LBAs. 

We also fit a log(L2-io kev) — lo gfLpTR) correlation from our LBA sample by applying 
the survival analysis (ASURV Rev 1.2; llsobe fc Feigelson Ill990l ). which gives 



log(L2_io kcv) = (0.41 ± 0.35) X log(LFiR) + (23.42 ± 11.28) 
The dash-dotted line shows the best fit for the correlation. 



(2) 



We find similar results in the soft X-ray band. We fit the correlation between the soft 
X -ray (0.5-2 keV ) and FIR luminosity of our LBA sample, and compared it with the result 
of iRanalli et al. I (120031 ). The survival analysis gave a free slope fit 



log(Lo.5-2 kev) = (0.42 ± 0.24) x log(LFiR) + (22.89 ± 10.76) 



(3) 



Besides the linear slope fit, iRanalli et al. I ( l2003l ) gave a free-slope fit of -Lo.5-2 keV oc 
-^FiR^° for star forming galaxies. In Fig. [HI the dash-dotted line indicates the best-fit of 
LBAs, and the Ranalli curve is shown as the dashed line below it. As in the case of the 
hard X-ray emission, only part of the soft X-ray fiux is likely to be contributed by normal 
processes associated with star formation. 

Given that the LBAs su ffer from less dust ext inction than typical local starbursts with 
similar star-formation rates ( lOverzier et al. 1 1201 if ), one might worry that the high ratio of 
X-ray to far-IR luminosity seen in Figures 7 and 8 is due to t he incomplete conversi on of 



far-UV into far-IR light. This is not the case. As discussed in lOverzier et al. I (120111 ). the 



star-formation rates derived from the observed far-UV are much smaller than those derived 
from the far-IR data, implying that the far-IR is a good proxy for the star-formation rate. 
Our results in section 3.2 and Figure 6 further support this. 



4.2. X-ray Properties compared to AGN 

The results from §4.11 are consistent with interpreting the optical spectra as indicating 
the presence of an AGN. The best tracer of both Type 1 and Type 2 AGN in the optical 
spectra is the [O III]A5007 emission-line. In this case, we would expect that the amount of 
hard X-ray emission from the LBAs would be consistent with the values found for AGN with 
similar [O III] luminosities. 

In Fig. ini we plot a histogram of the ratio of the hard X-ray (2-10 keV) luminosity 
obtained from the power-law model spectral fit of our LBAs to its observed [O III] luminosity. 
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We also show the empirical distribution of the lumi nosity ratios of Type 1 (dashed blue line) 



and Ty pe 2 AGN (dot-dashed red line) in this plot (IHeckman et al. Il2005bl : iLaMassa et al. 



2010bl ). The ratios of hard X-ray to [O III] luminosities in the LBAs are too low (by about 



an-order-of-magnitude) to be consistent with the values seen in Type 1 AGN. However, they 
do lie within the relatively broad range of values seen in Type 2 AGN. The relatively weak 
hard X-ray luminosity of Type 2 AGN is due to obscuration of the hard X-rays by the high 
column densities of gas associated with the AGN torus. In contrast, the [O III] line is formed 
in the 100 to 1000 pc-scale narrow line region and is not affected by the torus surrounding 
the black hole. 

At first sight, the presence of significant obscuration of the hard X-ray emission in the 
LBAs would not be consistent with the relatively low absorbing column densities listed in 
Table ID Here we see that the column density derived from the spectral fitting of an ab- 
sorbed power-law gives the 90% confidence level upper limit at the order of 10^^ cm~^ (10^^ 
cm~^ for SDSS J0802-I-3915). These columns are much too low to attenuate the hard X- 
ray emission by ~ an order-of-magnitude. Such a seeming paradox is often seen in Type 2 
AGN and is most likely due to fitting an overly idealized X-ray spectral mo del. More real 



istic r nodels (e.g. partial covering m odels and/or models with reprocessing) (ITurner et al. 



19971 : iLaMassa et al. 1 12009| . |2010b|) can reconcile the observed X-ray spectrum with heavy 



obscuration. 

Alternatively, the actual absorbing column densities derived from the fits to the X- 
ray spectra may be correct. In this case, the low ratio of hard X-ray to [O HI] luminosities 
relative to unobscured AGN would require that the [ O HI] line is pred o mina ntly produced by 
the starburst rather than the AGN. As discussed by lOverzier et al. I ( l2009l ). this possibility 
can not be excluded. 

Assuming for the moment that the relatively high ratios of X-ray to far-IR luminosity 
indicate that an AGN is present in these LBAs, our results still do not establish whether this 
AGN is more luminous and heavily obscured in hard X-rays or less luminous and relatively 
unobscured. An additional diagnostic of an obscured AGN is a strong iron Ka emi ssion- line 



with a large equivalent width of about one keV (e.g. ILaMassa et al. II2009I . l2010al ). In Fig. 
121 we show the comparison between the data and the spectral fit from an absorbed power- 
law. We find that there is a large data-to-model ratio at the Fe Ka energy in 1434+0207, 
0214-1-1259 and 0802+3915. We used a power-law plus Gaussian model to fit the PN spectra 
of 1434+0207 and 0214+1259 (insufficient photon counts for 0802+3915). We removed the 
5-8 keV photons first to determine the power-law index, then froze this parameter and added 
the 5-8 keV photons with a Gaussian component. The line energy is 6.6 keV for 1434+0207 
and 6.2 keV for 0214+1259, and the spectral plots are shown in Fig. [TOl The plots indicate 
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the possibility of an iron emission line in these LBAs. However, the photon counts are too 
few to give good estimates of the equivalent width of iron lines. 

In order to further investigate the presence of a Fe line in our sample, we loaded the 
spectra of all 6 LBAs in XSPEC to construct an average spectrum due to the small number 
of photons of each LEA (effectively stacking the spectra in model space). This was done in 
the same way by which we fitted the thermal component simultaneously in §3.1.21 This is 
not a fit to a stacked spectrum, but rather a joint fit to the set of observed spectra without 
shifting them to their rest-frame. The averaged spectrum was modeled as a continuum power- 
law in the range of 3-8 keV plus a Gaussian emission line. The intrinsic line width [a) in 
the Gaussian component was fixed at 0.01 keV, and the photon indices of the continuum 
power-law were set to 1.9. By tying the Gaussian normalization of all data sets, we obtain 
the emission line energy at 6.40 keV, with a 90% confidence region (AC-stat=4.6 for two 
model components) of 6.09-6.68 keV. The equivalent width (EW) is 1.0 keV from this joint 
fit, and the associated 90% confidence upper limit of EW is about 4.4 keV, where the upper 
limit is obtained by scaling the confidence regi on of the hne i ntensity by the best-fit value 



of the EW, which is likely to be overestimated ( lYaqoob 1119981 ). However, the data has low 
signal-to- noise ratio, and the emission line profile is not significant. When we removed the 
Gaussian component in the spectral fitting, the C-stat only increased by 5.5 compared to 
the value of 109.6/102(dof) for the model with Gaussian. 

We then performed the Markov Chain Monte Carlo (MCMC) simulation to test the 
significance of the iron emission line from the joint fit. We generated 6 chains, and the 
length of each chain is 5 x 10^. The 90% confidence level for the line energy from the MCMC 
is 6.10-6.93 keV, and the estimated upper limit of EW is 2.5 keV. The histogram distribution 
of the line energy and normalization parameters is shown in Fig. [TTl The significance of the 
Fe Ka line is calculated as the fraction of non-zero line normalization in the chains, which 
is found to be 95.3%, i.e., at the ~ 2a level. 



4.3. Mid-IR Constraints on a Type 2 AGN 



One other test for the presence of a Type 2 AGN in the LBAs is to compare the amount 
of mid-IR emission expe cted from a such an AGN with the observed [O HI]A5007 luminosity. 
LaMassa et al. I (j2010al ) have shown that both of these are good proxies for the bolometric 
luminosity of an obscured AGN (unlike the 2-10 keV X-ray luminosity). 

In Figure 12 we plot histogranas of the ratio of mid-IR/ [O HI] luminosities for the sample 
of Type 2 AGN in iLaMassa et al. I (l2010b( ) and for the LBAs in this paper. The distribution 



- 13 - 



for the Type 2 AGN has a mean of 1.89 dex and a dispersion of 0.57 dex compared to 2.71 
and 0.34 respectively for the LBAs. All six of the LBAs lie above the mean value for the 
Type 2 AGN. 

This suggests that if the LBAs are indeed starburst/AGN composite systems, the star- 
burst dominates over the AGN in the mid-IR in most cases. In a future paper we will 
examine this issue in more detail spectroscopically using Spitzer IRS data. Note that AGN 
emission from the torus peaks in the mid-IR, while the emission from dusty starbursts peaks 
in the far-IR (as is also the case for the LBAs). We therefore conclude that the bolometric 
luminosity of the LBAs is primarily due to the starburst. 

As noted above, it is possible that the [O III] emission in these LBAs is primarily 
produced by the starburst and that the hard X-ray source suffers little attenuation. In this 
case, the AGN contribution to the bolometric luminosity will be even smaller. 



4.4. Constraints on the Properties of Supermassive Black Holes 



As discussed in the introduction, these LBA-composites may offer us an unparalleled 
opportunity to study th e processes involved in the formation of sup ermassive black holes in 



the early universe (e.g., lElmegreen. Bournaud &: Elmegreen Il2008[ ). In this section we will 



infer the properties of the supermassive black holes in our sample if galaxies using standard 
methods. 

We will consider two possibilities. The first is the minimal case in which the hard X-ray 
so urce is not obscured. In this case, the bolometric correction to the hard X-ray luminosity 
in iMarconi et al I ( 120041 ) leads to a range of AGN bolometric luminosities from 1 to 6 x 10^ 
Lq. The implied black hole masses then range from ~ 3 x 10^ to ~ 2 x 10^ {L/LEdd)~^ M©. 
If these black holes are indeed radiating near the Eddington limit, the implied masses are 
then in the r ange of the "intermed iate mass" black holes found in some nearby low-mass 



galaxies (e.g. iPeterson et al. II2005I ). In the context of the high redshift universe, similar 



objects could represent the seeds for the subsequent growth of more massive systems. 

The second possibility is that the hard X-rays are significantly obscured in these LBAs 
and that the [O III]A5007 emission-line is powered by a Type 2 AGN. Adopting the bolo 
metric correction to the extinction-corrected [O III] luminosity from lKauffmann fc Heckman 
(I2OO9I ) then yields AGN bolometric luminosities of ~ 10^° to ~ 2 x 10^^ Lq and black hole 



masses of ~ 3 x 10^ to ~ 6 x 10^ {L/LEdd) ^ M©. In this case, an Eddington- limite d 



black hole would have a mass of-order 10 ^ of the DCO stellar mass (lOverzier et al 



This is similar to the mass ratio seen in present-day galactic bulges (IHaring &: Rix 



2009). 



2004 
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Marconi &: Hunt Il2003[ ). However, if the [O HI] line is contaminated or dominated by the 



starburst, the AGN bolometric luminosity and implied black hole mass would be overesti- 
mated. 



5. Conclusions 

We present in this paper XMM observations of 6 Lyman Break Analogs (LBAs): rare 
members of a local population that strongly resemble Lyman Break Galaxies at high-redshift 
in stellar mass, star formation rate, metallicity, size, dust extinction, and velocity dispersion. 
Our sample consists of LBAs whose optical spectra are consistent with a composite system 
consisting of both an intense starburst and a Type 2 (obscured) AGN. However, other in- 
terpretations of the optical spectra are possible. To further elucidate the nature of these 
objects we combine our new X-ray data with SDSS optical spectra and Spitzer mid- and 
far-IR photometry. We summarize the key results as follows: 

(1) We fitted their X-ray spectra with an absorbed power- law model with the parameters 
typically found in AGN. The column densities found from the spectral fit are less than 
1022 cm-2. 

(2) We found that the ratios of X-ray to FIR luminosity of the LBAs are higher than the 
typical starburst values by factors of three to thirty. This implies a source of X-ray emission 
in addition to what is seen in typical starbursts. We suggest that the most likely candidate 
would be an AGN. 

(3) The ratios of hard X-ray (2-10 keV) to [O HI] luminosity are about an order-of- 
magnitude smaller than in unobscured (Type 1) AGN but consistent with the broad range 
of values seen in Type 2 AGN. This might suggest that the hard X-rays suffer a significant 
amount of attenuation (as in the case of typical Type 2 AGN). This would imply that the 
fitted column densities derived from the simple absorbed power-law model underestimate the 
true attenuation (as is often the case in Type 2 AGN). An alternative interpretation is that 
the hard X-rays are relatively unobscured and that the [O HI] line is primarily produced by 
an extreme starburst. 

(4) To discriminate between these possibilities (and provide a confirmation of the pres- 
ence of an obscured AGN), we searched for the presence of the 6.4 keV Fe Ka emission line 
in our LBA sample. However, the limited photon detections and low signal-to-noise ratio do 
not allow us to come to a strong conclusion. The MCMC simulation gives a ~ 2cr significance 
level for the iron line and a 90% confidence upper limit of 2.5 keV for the equivalent width. 



- 15 - 



(5) We find that the ratio of mid-IR (24/im) and [O III]A5007 luminosities of the LBAs 
are higher than the values seen in Type 2 AGN (by a difference of 0.8 dex in the mean). 
This would imply that the bolometric luminosity of these systems is primarily due to the 
starburst rather than an AGN. If the [O III] emission is primarily due to a starburst, this 
conclusion is even stronger. 

(6) Given that these objects have significantly higher ratios of X-ray to far-IR luminosity 
than starbursts, we conclude that the most likely interpretation is that these galaxies are 
indeed composites of an intense starburst and an AGN (consistent with their optical emission- 
line spectra). If so, the implied black hole masses are of-order 10^ for Eddington- 
limited accretion in the case where the hard X-rays are unobscured and roughly an order- 
of-magnitude larger for an obscured AGN that produces the observed [OIII] emission-line. 
These may be good analogs to black holes forming in the dense, stellar clumps in galaxies 
in the early universe. 

We are very grateful to the anonymous referee for helpful comments that have improved 
the manuscript. We also thank Marat Gilfanov for discussion. This work is supported by 
NASA grant NNX08AZ0G. 
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Table 1: XMM- Newton observations of Lyman Break Analog targets 



Name 


Observation ID 


Date 


Exp. Time^ A^hg'' 


z 


Mode 


Filter 


(SDSS J...) 




(DD/MM/YY) 


(ks) 








080844.26+394852.3 


0553790101 


03/10/08 


8.7/11.5/11.5 5.05 


0.091 


FF^ 


medium 


210358.74-072802.4 


0553790201 


18/05/08 


8.7/11.5/11.5 5.41 


0.137 


FF 


thin 


143417.15+020742.3 


0553790301 


20/07/08 


13.7/22.1/22.1 2.80 


0.180 


FF 


thin 


021348.53+125951.4 


0553790401 


08/01/08 


19.4/23.6/23.7 7.06 0.219 


FF 


thin 


092159.38+450912.3 


0553790501 


27/04/09 


3.7/20.6/22.4 1.55 


0.235 


FF 


medium 


080232.34+391552.6 


0553790701 


24/03/09 


1.3/6.4/6.3 5.15 


0.267 


FF 


thin 



^The exposure time is the good time interval after flare filtering. We list the exposure time for all three 
detectors (PN/M0S1/M0S2). 

''Galactic column densities from the GOLDEN tool,in units of 10^*^ cm~^. 
^FuU frame imaging window mode. 
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Table 2. Ultra- Violet observations of Lyman Break Analog targets. 





TD 


Piltpr 

J. 11 LCI 


V^VJllC^^tCLl IVcttC 


A/r f» crn 1 n H p 

iVlct^lll t LILIC 


UV flii-x- 

VJ V 11 Ll-A. 


Luvwi 








s-1 




(10-^6 erg s-^ cm-2 A"^) 


lO^^^erg s-^ 


SDSS 


J0808+3948 


UVWl 


2.220 ±0.092 


16.34 ±0.045 


10.57 ±0.44 


2.70 


SDSS 


J2103-0728 


UVWl 


0.938 ±0.038 


17.28 ±0.045 


4.45 ±0.18 


2.74 


SDSS 


J1434+0207 


UVWl 


0.475 ±0.039 


18.01 ±0.089 


2.26 ±0.19 


2.09 






UVM2 


0.141 ±0.024 


17.90 ±0.188 


3.10 ±0.53 




SDSS 


J0214+1259 


UVWl 


0.502 ± 0.054 


17.95 ±0.116 


2.39 ±0.25 


3.85 






UVM2 


0.099 ± 0.027 


18.28 ±0.291 


2.18 ±0.59 




SDSS 


J0922+4509 


UVWl 


0.882 ± 0.045 


17.34 ± 0.056 


4.20 ±0.21 


7.93 






UVM2 


0.265 ±0.025 


17.21 ±0.104 


5.83 ±0.55 








UVW2 


0.119 ±0.020 


17.18 ±0.178 


6.79 ± 1.14 




SDSS 


J0802+3915 


U 


0.587 ±0.033 


18.84 ±0.061 


1.14 ±0.06 








UVWl 


0.252 ±0.033 


18.84 ±0.061 


1.20 ±0.16 


3.03 






UVM2 


0.114 ±0.019 


18.13 ±0.183 


2.51 ±0.42 








UVW2 


0.056 ±0.018 


17.99 ± 0.338 


3.20 ±1.03 





Note. — The filters with decreasing wavelength are U (344nm), UVWl (291nm), UVM2 (231nm) 
and UVW2 (212nm), and the corresponding count rates/fiux conversion factors are based in white dwarf 
standard stars. 



Table 3: [O III] and IR luminosities 



ID 


-^[OIII],obs 


-^[OIII],corr 




-^FIR 




(10^1) 


(10^1) 


(10^^) 


(10^4) 


SDSS J0808+3948 


0.43 ±0.01 


0.66 


0.76 


1.1 ±0.2 


SDSS J2103-0728 


5.00 ±0.06 


12.03 


4.8 


8.8 ± 1.3 


SDSS J1434+0207 


1.32 ±0.03 


2.84 


0.48 


2.0 ±0.6 


SDSS J0214+1259 


0.57 ±0.04 


1.22 


1.9 


9.1 ±1.7 


SDSS J0922+4509 


2.05 ±0.06 


4.75 


2.4 


14.3 ±2.2 


SDSS J0802+3915 


1.54 ±0.05 


4.47 


1.2 


6.5 ±2.0 



Note. — [O III] luminosities are obtained from MPA/JHU cata log of SDSS DR7. The mid- and far-IR 
luminosities are based on Spitzer photometry (jOverzier et al. 11201 if ). All luminosities are in units of erg s~^. 



Table 4. X-ray spectral fits using the C-statistic. The errors and upper limit correspond to the 90 per cent 

confidence level. 



ID Counts 


r C-stat/dof A^H.in (r = 


1.9) C-stat/dof 


.^^0.5-2keV* 


-^2-lOkeV* 


.^^2-lOkev/ -^^ [GUI], corr 


PN/M1/M2 




(lO^^cm- 




(10" erg s-i) 


(10^1 erg s-i) 




(1) (2) 


(3) (4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SDSS J0808+3948 75/34/44 


1.66 (1.41-1.92) 116.11/138 


< 0.14 


115.82/138 


2.11 (1.55-2.80) 


4.01 (2.94-5.33) 


6.08 


SDSS J2103-0728 76/26/17 


2.48 (2.08-2.94) 112.78/100 


< 0.04 


118.60/100 


2.93 (2.14-3.95) 


4.83 (3.53-6.51) 


0.40 


SDSS J1434+0207 86/41/43 


1.46 (0.88-2.04) 134.41/156 


< 0.27 


135.80/156 


2.65 (1.53-4.28) 


4.41 (2.55-7.13) 


1.55 1 



SDSS J0214-M259 83/25/28 2.72 (1.91-3.68) 133.53/122 < 0.09 136.31/122 2.39 (1.48-3.68) 4.30 (2.66-6.62) 3.51 to 

SDSS J0922-I-4509 54/70/62 1.89 (1.50-2.30) 172.57/165 < 0.23 171.02/165 12.35 (6.64-19.94)22.70 (12.21-36.66) 4.78 | 

SDSS J0802-I-3915 15/8/19 1.71 (-0.58-4.17) 28.56/38 < 1.37 28.62/38 9.29 (0.81-31.96) 16.80 (1.47-57.81) 3.76 



Note. — Column (1): LBA names; column (2): numbers of counts detected by PN, MOSl and M0S2; column (3): derived photon indices (F), 
no intrinsic absorption assumed; column (4) and (6): reduced C-statistic; column (5): derived intrinsic column densities at the redshifts of the 
sources: by freezing F at 1.9; column (7): luminosities in soft X-ray band (0.5-2.0 keV); column (8): luminosities in hard X-ray band (2.0-10.0 

keV); column (9): ratio of hard X-ray and extinction-corrected [O 111] luminosities. 

*The errors of luminosities correspond to the errors of normalizations of power-law components in the spectral fit. 
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Table 5. X-ray spectral fits using the x^-statistic. The errors and upper hmit correspond 

to the 90 per cent confidence level. 



ID 


r 




7VH,i, (for r = 1.9) 










(1022cm-2) 




(1) 


(2) 


(3) 


(4) 


(5) 


SDSS J0808+3948 


1.71 (1.33-2.13) 


1.01 


< 0.24 


0.553 


SDSS J2103-0728 


3.19 (2.40-4.28) 


0.19 


< 0.08 


2.15 


SDSS J1434+0207 


1.24 (0.40-2.06) 


1.09 


0.16 (0-2.44) 


1.17 


SDSS J0214+1259 


> 2.05 


0.79 


< 0.68 


1.22 


SDSS J0922+4509 


2.83 (2.01-4.18) 


0.72 


< 0.18 


0.99 



Note. — Column (1): LBA name; column (2): derived photon indices 
(r), no intrinsic absorption assumed; column (3) and (5): reduced x^; col- 
umn (4): derived intrinsic column densities at the redshifts of the sources; 
by freezing F at 1.9. 



Table 6. X-ray spectral fits with both photon indices and column densities free 
parameters. The errors and upper limit correspond to the 90 per cent confidence level. 



ID 


statistic 


r 


A^H,in 


xl or 




method 




(1022cm-2) 


C-stat/dof 


SDSS J0808+3948 




1.85 (1.22-3.05) 


< 0.53 


0.69 




c 


1.79 (1.35-2.50) 


< 0.22 


0.86 


SDSS J1434+0207 


x' 


1.23 (0.02-3.28) 


< 21.8 


1.18 




c 


1.41 (0.713-2.47) 


< 0.35 


0.87 


SDSS J0922+4509 


x' 


4.56 (1.87-10.0) 


0.59 (0-2.91) 


0.71 




c 


3.35 (1.94-6.56) 


0.42 (0.05-1.41) 


1.01 



Note. — The errors of photon indices and column densities correspond to Ax^ 
or AC=4.6, since there are two interesting parameters. 



Table 7. apec plus power-law model 



ID 


Temperature 
(eV) 


ATh 
(1022cm-2) 


-^thermal 

(10-*i erg s-i) 


T 2 
^0.5-2keV,PL 

(10*1 erg s-i) 


i2-10keV"^ 

(10"! erg s-i) 


C-stat/dof 


SDSS J0808+3948 


500 (fixed) 


0.10 (0.01-0.41) 


0.29 (0.15-12.69) 


1.78 (1.09-2.74) 


3.76 (2.31-5.79) 


114.8/133 


SDSS J2103-0728 


422 (276-755) 


0.17 (0-1.03) 


1.74 (0.16-3.50) 


1.45 (0.56-2.63) 


3.85 (1.49-6.75) 


101.9/96 


SDSS J1434+0207 


339 (217-689) 


1.07 (0.18-2.55) 


1.48 (0.39-2.98) 


0.86 (0.19-1.75) 


6.61 (1.47-13.44) 


129.3/152 


SDSS J0214+1259 


570 (326-844) 


266.5 (13.9-) 


2.23 (1.05-3.84) 


1.95 (1.0-4.73) 


6.20 (3.10-15.03) 


125.7/118 


SDSS J0922+4509 


500 (fixed) 


0.16 (0-0.65) 


2.98 (1.50-11.10) 


6.04 (1.13-12.51) 


21.84 (4.10-45.23) 


169.4/162 


SDSS J0802+3915 


500 (fixed) 


1.96 (0.81-50.35) 


12.64 (1.97-31.67) 


14.23 (0.23-166.63) 


40.58 (0.66-475.17) 


24.42/35 



^Thermal emission luminosity. 
2power-law luminosity in 0.5-2 keV band. 
^Total luminosity in 2-10 keV. 



Table 8. apec plus power-law model in a joint fit 



ID A^H -f'thcrmal^ -f'0.5-2keV,PL^ -f'2-lOkeV^ 

(1022cm-2) (10^^ erg s"^) (10^^ erg s"^) (10^^ erg s"- 



SDSS J0808+3948 
SDSS J2103-0728 

SDSS J1434+0207 
SDSS J0214+1259 
SDSS J0922+4509 
SDSS J0802+3915 



0.14 (0.07-0.28) 
0.35 (0.05-1.17) 

0.59 (0.06-1.86) 
769 (45-1240) 
0.19 (0.05-0.52) 
2.94 (0.50-9.30) 



0.39 (0.13-0.75) 
1.83 (0.60-3.48) 
1.12 (0.37-2.12) 
1.94 (0.64-3.69) 
1.27 (0.42-2.42) 
6.34 (2.08-12.05) 



2.18 (1.33-3.36) 
0.97 (0.37-1.76) 
1.57 (0.35-3.19) 

8.29 (1.55-17.17) 
2.31 (0.04-27.05) 



4.87 (2.98-7.50) 
3.31 (1.28-6.00) 
6.91 (1.53-14.06) 
9.11 (4.67-22.1) 
19.43 (3.64-40.24) 
63.49 (1.03-743.45) 



^Thermal emission luminosity. 
^Power-law luminosity in 0.5-2 keV band. 
^Total luminosity in 2-10 keV. 
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Table 9. SFRs derived from UV luminosities. 



ID 


SFRfUVWl) 

kJ J- J- ^ V ^ V V V -L J 


SFRfH -\-2Aam)* 




Moyr-i 


Moyr-i 


SDSS J0808+3948 


6.5 


16.1 


SDSS J2103-0728 


6.7 


108.3 


SDSS J1434+0207 


5.1 


20.0 


SDSS J0214+1259 


9.3 


35.1 


SDSS J0922+4509 


26.1 


55.1 


SDSS J0802+3915 


7.3 


30.4 



*The values of SFRs derived from combination of the 
uncorrected Hp, luminosity and 2Af xm luminosity are ob- 



tained from lOverzier et al. I (120091 ) 
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Fig. 1. — Optical emission line diagnostic d iagram showing the 31 w e ll-studied nearby (z < 



0.3) LB As (colored symbols ) discovered by lHeckman et al. I (l2005al ): iHoopes et al. I ( 120071 ). 
The LBAs marked in red in this paper lie in the star-forming/AGN composite region (we 
did not get obseryation time of the object marked i n green), where the s olid line is from 



Kauffmann et al. I (l2003l ) , and the dashed line is from iKewley et al. I (l2006l ) 
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Fig. 2. — X-ray spectral plots of each LBA-composite fitted with C-statistic, and the ratios 
between the model and data of the fit are shown in the bottom panels. The colors of black, 
red and green indicate the data in detectors of PN, MOSl and M0S2, respectively. 
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Fig. 3. — X-ray spectral plots of each LBA-composite fitted with x^-statistic, and the ratios 
between the model and data of the fit are shown in the bottom panels. The colors of black, 
red and green indicate the data in detectors of PN, MOSl and M0S2, respectively. The plot 
for SDSS J0802-I-3915 is not shown due to insufficient counts for x^-statistic. 
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Fig. 4. — Comparison of contributions to the soft X-ray luminosity between power-law and 
thermal emission components in the power-law plus thermal spectral fits. Left panel shows 
the results from the individual fits (Table [7]), and the right panel shows those from the joint 
fit (Table |8]). The error bar for soft X-ray luminosity is not shown as it makes the plot 
unwieldy. SDSS J0214+1259 has a negligible 0.5-2 keV power-law luminosity due to its 
Compton-thick obscuration estimated in the simultaneous joint fit, so it is not included in 
the right panel. The dashed line indicates where both values equate. 



- 29 - 




Fig. 5. — UV images (UVWl filter) of eacli LBA-composite target. Tlie size of each image 
is 0.6'x0.6'. 
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Fig. 6. — Star formion rates derived from the UV luminosity and from Rq. + 24/xm emission 
lines. The units are per year. 
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Fig. 7. — Hard X-ray luminosity versus far infrared luminosity. Our sample of LBA- 



composites is shown as plus sign s with error bars. Others are from llwasawa et al. I (120091 ) 



(Iw09) and iLehmer et al. I (120101 ) (LehmerlO) for comparison, where the diamond symbols 
indicate the Iw09 AGN samples, and asterisks and upper limit arrows are "Hard X-ray 
Quiet" galaxies (HXQ) from Iw09, and LIRG samples from LehmerlO are marked as tri- 
angles (triangles in red are identified as AGN). The solid line indicates th e empirical hard 
X-ray and far infrared luminosity correlation given by lRanalli et al. I (120031 ) . and the dashed 
lines show an associated uncertainty of a factor of 2. The dash-dotted line is the correlation 
derived from our LBA sample by using the survival analysis. 
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Fig. 8. — Soft X-ray (0.5-2.0 keV) lu minosity versus far infr ared luminosity. The dashed line 
is the -Lo.5-2.okcV — -^^fir relation from lRanalli et al. I ( 120031 ) . while our best-fit correlation is 
indicated by the dash-dotted line. 
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Fig. 9. — Histogram of the log of -^2-10keV/-^[OIII] obs ^^^^^s- sample is represented 
by the solid black li ne. The dashed b l ue line and dot-dashed red line show the typical Type 
1 and 2 AGNs from lfleckman et al. I fconSbl ). 
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Fig. 10. — Evidence for iron emission lines on power-law continuum spectra for 
SDSS J1434+0207 and SDSS J0214+1259. 
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Fig. 11. — Histogram distribution of the line energy and normalization parameters from the 
MCMC simulation. 
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Fig. 12. — Histograms of the log of the ratio of the mid- IR (24Atm) and extinct ion-corrected 
[O III]A5007 emission-hne luminosity for Type 2 AGN (ILaMassa et al. Il2010bh (black solid 
line) and the current sample of LBA-composites (red dashed line). The ratio is systematically 
higher in the LBA-composites (2.71±0.34) than that in Type 2 AGN (1.89±0.57), suggesting 
that the mid-IR continuum is primarily due to the starburst in most of the LBA-composites. 



